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Expansion of the Homeostasis Model 
Assessment of p-Cell Function and 
Insulin Resistance to Enable Clinical 
Trial Outcome Modeling Through the 
Interactive Adjustment of Physiology 
and Treatment Effects: IHOMA2 



OBJECTIVE — To describe and make available an interactive, 24-variable homeostasis model 
assessment (iHOMA2) that extends the HOMA2 model, enabling the modeling of physiology 
and treatment effects, to present equations of the HOMA2 and iHOMA2 models, and to exem- 
plify iHOMA2 in two widely differing scenarios: changes in insulin sensitivity with thiazolidi- 
nediones and changes in renal threshold with sodium glucose transporter 2 (SGLT2) inhibition. 

RESEARCH DESIGN AND METHODS— 1HOMA2 enables a user of the available soft- 
ware to examine and modify the mathematical functions describing the organs and tissues in- 
volved in the glucose and hormonal compartments. We exemplify this with SGLT2 inhibition 
modeling (by changing the renal threshold parameters) using published data of renal effect, 
showing that the modeled effect is concordant with the effects on fasting glucose from indepen- 
dent data. 

RESULTS — iHOMA2 modeling of thiazolidinediones effect suggested that changes in insulin 
sensitivity in the fasting state are predominantly hepatic. SGLT2 inhibition modeled by iHOMA2 
resulted in a decrease in mean glucose of 1.1 mmol/L. Observed data showed a decrease in 
glucose of 0.9 mmol/L. There was no significant difference between the model and the indepen- 
dent data. Manipulation of iHOMA2's renal excretion threshold variable suggested that a de- 
crease of 17% was required to obtain a 0.9 mmol/L decrease in mean glucose. 

CONCLUSIONS — iHOMA2 is an extended mathematical model for the assessment of insulin 
resistance and P-cell function. The model can be used to evaluate therapeutic agents and predict 
effects on fasting glucose and insulin and on p-cell function and insulin sensitivity. 
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Type 2 diabetes is caused by a combi- 
nation of progressive P-cell dysfunc- 
tion, relative insulin deficiency, and 
variable degrees of insulin resistance that 
lead to dysregulation of glucose homeo- 
stasis. Understanding the biochemistry, 
phenotypic details, and genetic mecha- 
nisms contributing to this can yield im- 
portant information on pathophysiology. 
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The progressive nature of the disease, as 
well as measuring the rate of deteriora- 
tion, has presented an ongoing challenge 
to clinicians and scientists alike. Tools to 
track P-cell functional changes and insu- 
lin resistance fall into three broad catego- 
ries: measures of glycemic status (e.g., 
fasting glucose, HbA lc ), physiological in- 
vestigations (e.g., clamp techniques [1,2], 



glucose tolerance tests), and mathemati- 
cal modeling (e.g., minimal model [3], 
Mari model [4,5], homeostasis model as- 
sessment [HOMA] [6-9]). 

No single approach proved sufficient, 
either, for a comprehensive quantitative 
description of P-cell dysfunction or 
insulin resistance. Measures of these pa- 
rameters vary depending on whether 
measurements are from basal or stimu- 
lated or fasting or postprandial subjects 
and whether pharmaceutical agents are 
being taken. Physiological techniques, 
ranging from simple glucose tolerance tests 
to euglycemic clamps and stable isotope 
studies, require expertise and are time 
and resource intensive, limiting their use 
to relatively small numbers of subjects (10). 

Mathematical modeling techniques 
also vary in their physiological assump- 
tions. Computer-based solutions from 
clinical interventions (e.g., oral glucose 
tolerance tests with "minimal model" 
readout) have limitations because of the 
high number of samples required from 
each subject. Simpler modeling methods 
(e.g., HOMA2) use paired fasting plasma 
insulin and glucose concentrations to de- 
rive data on P-cell function and insulin 
sensitivity. HOMA2 yields a single read- 
out of p-cell function and insulin resis- 
tance for each subject and has the 
advantage that, since it only requires 
paired basal insulin and glucose measure- 
ments, it can be used in large epidemio- 
logical and pharmaceutical studies. 

One disadvantage with HOMA2 is 
that it is not an appropriate model to use 
when evaluating treatments that have 
similar functional effects on blood glu- 
cose but different modes of action. For 
example, in HOMA2 P-cell function is 
characterized internal to the model as a 
sigmoidal dose response curve relating 
insulin secretion to the prevailing glu- 
cose concentration. The shape of this 
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sigmoidal curve is modeled using two 
principal variables, one of which de- 
scribes the rate of insulin secretion (K m ) 
and one that describes the maximal insu- 
lin secretion (V max )- In HOMA2, changes 
in p-cell function, altering insulin secre- 
tion, are entirely attributed to changes in 
V max and K m in the basal state (11). This 
modeling of (3-cell function has been a 
good approximation, but a comprehen- 
sive description of the p-cell dynamics re- 
quires more variables (11). The (3-cell 
functional changes caused by sulfonyl- 
ureas (12) have different dose response 
changes and characteristics from those in- 
duced by incretin hormones (13). Similarly, 
total body insulin sensitivity, defined as net 
glucose clearance for any given insulin 
concentration, cannot be comprehen- 
sively described using one fixed function 
as it is in HOMA2. Hepatic insulin sensi- 
tivity may differ from peripheral sensitiv- 
ity, and change in glucose clearance can 
be completely unrelated to either hepatic 
or peripheral sensitivity as in the case of 
sodium glucose transporter 2 (SGLT2) 
inhibition (14). The invariant nature of 
HOMA2 does not allow a sufficient descrip- 
tion for other organs and tissues involved in 
glucose homeostasis, and this can be impor- 
tant where there is knowledge of changed 
function, e.g., in the liver, which could be 
used to improve model outputs. 

We have developed an interactive, 
24-variable model (termed iHOMA2) that 
addresses problems associated with the 



fixed assumptions within HOMA2. Here, 
we describe the model in cartoon, graphical, 
and mathematical detail; discuss how it can 
be used and manipulated; present an ex- 
plicit example of its use to model changes 
in renal threshold induced by blockade of 
glucose reabsorption by SGLT2 inhibition; 
and show how this affects the modeled 
fasting glucose. The model is available on 
open access. 

RESEARCH DESIGN AND 
METHODS 

Model development 

The iHOMA2 model is shown in Fig. 1 as 
graphical (A), box diagrammatic (B), and 
mathematical (C), respectively. iHOMA2 
is an integrated computer-based mathe- 
matical model of glucose and hormonal in- 
teraction under homeostatic conditions. 
The model, now available online at http:// 
www.ihoma.co.uk, runs in real time with 
24 operator-controlled variables (Table 1) 
and graphical output displays. The baseline 
characteristics were built from those used 
in the original H0MA2 model, with all of 
the dose-response variables now explicit. 
iH0MA2 runs interactively and exactly 
for each calculation. iH0MA2 in its default 
start-up setting gives identical readings to 
H0MA2 and can be used as a direct sub- 
stitute for H0MA2 in this mode. The op- 
erator can modify each of the variables 
using an interactive sliding control display. 
The operator can control every aspect of the 



dose-response curve. For example, the 
(3-cell characteristics are described by 
P1-P5, each of these being independently 
adjustable. This allows "what if scenarios 
to be explored: "What would be the effect 
on glucose if V max of p-cell function were 
50%?" "How might that be modified if the 
dose response curve were shifted to the left?" 
"What if autonomous insulin secretion 
continued at low blood glucose?" Simi- 
larly, the functions relating to the other 
organs and tissues involved in the glucose 
and hormonal compartments can be mod- 
ified using sliding control displays. In the 
H0MA2 model, insulin sensitivity is treated 
as a whole-body effect, altering the liver 
and periphery to the same extent. In 
iH0MA2, this has been uncoupled and 
the insulin sensitivity of these organs and 
tissues can be modified independently. The 
ability to alter the 24 variables of iH0MA2 
enables the modeling of known or surmised 
pathology and physiology and the effect of 
treatments both alone and in combination. 
The effects of the treatments on fasting glu- 
cose, insulin, (3-cell function (%B), and in- 
sulin sensitivity (%S) are graphically 
represented in the model. 

The model allows for analytical and 
predictive modes of use. The analytical 
mode allows insulin resistance and (3-cell 
function to be read from the input of in- 
sulin and glucose in the basal state, while 
the predictive function shows the esti- 
mated and modeled insulin and glucose 
concentrations in the basal state when the 
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eg = g * (PI A bkexp), bkexp = 0.25 , 

u = peripheral insulin concentration * PE5, 

uh = hepatic insulin concentration * L5, kGrad = Kidney gradient of slope, 
Kint = Kidney intercept of slope 



Figure 1 — The iHOMA2 model is represented in three different ways: A) as a cartoon displaying the organs and tissues involved; B) showing the dose 
response curves as graphs modeled for each organ or tissue; and C) providing the equations that govern the organs and tissues functioning in the 
model. B 2 and B 2 are modifiable for the brain, PE 1 -PE 5 /or the periphery , Rj-R 4 for the kidneys, G 1 and G 2 for the gut, Lj-L 5 for the liver, and Pj— P5 
for the pancreas (see Table 1). 
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iH0MA2: an interactive HOMA model 
Table 1 — Modifiable variables within iH0MA2 





Name 


Description 


Pancreas 


Pi 




|3-Cell function of the islets 


Pi 


Max 


Maximum output of insulin 


P 3 


Effective concentration 


Glucose concentration required 


?4 


Slope factor (Hill slope) 


Rate of insulin secretion modifier 1 




Liver 


L 2 


Insulin effect 


Insulin concentration effect 


i-3 


Max 


Maximum glucose output 


u 


Slope 


Rate of glucose output 


u 


%S 


Sensitivity of glucose output to insulin 


Periphery (muscle, fat) 


PEi 


Max 


Maximum glucose uptake 


PE 2 


Slope A 


Primary rate of glucose uptake 


PE 3 


Slope B 


Secondary rate of glucose uptake 


PE 4 


Insulin effect 


Insulin concentration effect 


Gut 


Gi 


Max 


Maximum glucose uptake 


G 2 


Slope 


Rate of glucose uptake 




Bi 


Max 


Maximum glucose uptake 




Renal 




R 2 


Slope B 


Secondary rate of glucose excretion 


R3 


Delay 


Glucose concentration required before 
glucose excretion started 


R 4 


Threshold 


Glucose concentration basal threshold 



The iHOMA2 variables that can be modified by the operator either by using a sliding scale on the interface or 
by entering exact values. The absolute values at default are equal to those used in HOMA2. 



p-cell function and insulin resistance pa- 
rameters are set. 

This article shows two detailed quan- 
tified scenarios to illustrate the interactive 
modalities. The first example shows that 
the effect of pioglitazones (thiazolidine- 
diones) on insulin resistance can be par- 
titioned between the liver and periphery. 
The second example illustrates the mod- 
el's use elucidating the effect of an SGTT2 
inhibitor on glycemia. All analyses were per- 
formed using SPSS, version 19.0 (SPSS, Chi- 
cago, IT). Statistical comparisons were made 
using Z tests for skewness, Student indepen- 
dent samples t test for comparison of means, 
and F tests for assessment of fit of the model 
to the observed data (15). 

Quantitative model usage: effect of 
pioglitazone 

To model the effects of pioglitazone, we 
examined the outcome when insulin sen- 
sitivity was modeled to be in the liver, in 



the periphery, or at both sites equally — all 
with a standardized increase in (3-cell 
function. Three possible sites of action 
on insulin sensitivity for pioglitazone 
were modeled: hypothesis 1, insulin sen- 
sitivity increases in both periphery and 
hepatic (variables L 5 and PE 5 ); hypothesis 
2, hepatic insulin sensitivity increases and 
peripheral insulin sensitivity remains un- 
changed (variable L 5 ); and hypothesis 3, 
hepatic insulin sensitivity remains un- 
changed and peripheral insulin sensitivity 
increases (variable PF 5 ). 

The development dataset comprised 
insulin and glucose values from a mono- 
therapy study of pioglitazone (16). The 
changes to (3-cell function and insulin 
sensitivity observed between the baseline 
visit and end of study (12 months later) in 
the development group were used as in- 
puts to adjust the variables in the 
iHOMA2 model for each of the three hy- 
potheses. A separate study of pioglitazone 



(17) was used for the verification group. 
The data from the baseline visits of the 
verification group were submitted to the 
adjusted iHOMA2, using the model in 
predictive mode, to determine the effect 
of pioglitazone for each of the three hy- 
potheses using as output the expected 
fasting glucose and insulin after therapy. 
We assessed bias and agreement using a 
Bland-Altman plot and assessed the fit of 
the model by examination of the least 
squares deviation from the line of unity 
(where the observed values equal the pre- 
dicted values) using an F statistic to test 
the model fit. 

Quantitative model usage: the effect 
of glucose reabsorption inhibitors 

SGTT2 partially prevent the reabsorption 
of glucose, thus changing the renal 
threshold for glycosuria. Patients using 
an SGTT2 inhibitor demonstrate a 
marked increase in glucose urinary loss. 
A recent publication estimated, for a 5-mg 
dose of dapagliflozin over 2 weeks, a 20% 
decrease in the renal threshold for gly- 
cosuria (18). To model the effect of 
SGTT2 inhibition, we used a phenotypi- 
cally similar subject set (previously pub- 
lished [19]) where fasting glucose and 
insulin measurements were known. SGLT2 
inhibition was modeled in iHOMA2 (vari- 
able R4 = 120%) based on the published 
data to predict the effect on the fasting 
glucose and insulin. The glucose changes 
predicted were compared with a separate 
preselected phase 3 trial of dapagliflozin 
(20) using an independent t test. 

Further, the iHOMA2 model was 
used to examine any change required in 
renal glucose excretion to achieve an 
equivalent change in fasting glucose in 
our cohort. iHOMA2 was manipulated 
(by change of variable R 4 ) until the change 
in glucose from the model was equal to 
the change in glucose from the dapagliflozin 
study. The value of the variable £4 was then 
read from the model. 

RESULTS 

Quantitative model usage: effect of 
pioglitazone 

Data were analyzed for normal distribu- 
tion before parametric analysis. The z test 
value for the skewness of the insulin data 
was 3.9. Tog transformation eliminated 
the skewness as indicated by the z value 
of 0.5. Geometric mean insulin results are 
therefore presented. 

In the development group, pioglita- 
zone increased p-cell function from 36.9 
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Table 2 — Example of hypothesis predictions from the iHOMAl model 



Insulin sensitivity 


Mean Aglucose 
(mmol/L) 




Geometric mean Ainsulin 
(pmol/L) 


Observed Modeled 


F 


Observed Modeled F 


Hypothesis 1: 50% hepatic, 








50% peripheral 


-2.4 -2.7 


143* 


-8.2 -14.6 970* 


Hypothesis 2: 100% hepatic 


-2.4 -2.6 


139* 


-8.2 -11.5 259* 


Hypothesis 3: 100% peripheral 


-2.4 -2.9 


141* 


-8.2 21.4 650* 



Example 1: Modeling drug effects (pioglitazone). Change in mean glucose and insulin predicted by the 
iHOMA2 model compared with observed results from an independent study using the F statistic to test model 
fit. *Signincant to P < 0.001. 



iHOMAl: an interactive HOMA model 

to 49.2, a relative increase of 33.4%, and 
increased insulin sensitivity from 61.4 to 
79.3, a relative increase of 29.1%. Since 
we now had a calculated generic change 
in insulin sensitivity with an associated 
p-cell functional change, we were able 
to partition the site of resistance to test 
three hypotheses (Fig. 2): 

1. That insulin resistance change with 
pioglitazone is equally partitioned be- 
tween the liver and the periphery. 
Whence we set hepatic insulin sensi- 
tivity increase by 29.1% and peripheral 
insulin sensitivity increase by 29.1% 
with |3-cell function increased by 33.4% 
(variables L 5 = 129.1%, PE 5 = 129.1%, 
P± = 133.4%). 

2. That insulin resistance change with 
pioglitazone is only sited at the liver. 
Whence we set hepatic insulin sensi- 
tivity increase by 58.2%, peripheral 
insulin sensitivity unchanged, and 
P-cell function increase by 33.4% 
(variables L 5 = 158.2%, PE 5 = 100.0%, 
P 1 = 133.4%). 

3. That insulin resistance change with 
pioglitazone is only sited at the periph- 
ery. Whence we set peripheral insulin 
sensitivity increase by 58.2%, hepatic 
insulin sensitivity unchanged, and 
p-cell function increase by 33.4% (var- 
iables L 5 = 100.0%, PE 5 = 158.2%, P 1 = 
133.4%). 

The predicted change in glucose and 
the change in insulin from these hypoth- 
eses were contrasted with the observed 
values from the verification group using 
Bland-Altman plots and F tests for model 
fit. The Bland-Altman plots did not show 
significant differentiation between the 
models. F test (Table 2) demonstrated 
that the second hypothesis — maximal he- 
patic insulin sensitivity — was found to be 
the model with the best fit to the observed 
data. The F statistics for linearity were 139 
and 259 (with 668 df) (significance P < 
0.001) for glucose and insulin, respec- 
tively. 

Quantitative model usage: the effect 
of glucose reabsorption inhibitors 

iHOMA2 predicted that a decrease by 
20% in the renal excretion threshold 
using 5 mg dapagliflozin treatment over 
2 weeks would result in a decrease in 
mean (SD) glucose of 1. 1 (0.3) mmol/L in 
the study subjects. Data extracted from 
published literature on dapagliflozin (5-mg 
dose) reported a similar change in mean 
glucose of 0.9 mmol/L in their subjects after 



2 weeks. There was no significant differ- 
ence between the model and the observed 
data. There are no published data on 
changes in insulin, so these could not be 
compared. 

Although the model and observed 
results were broadly concordant, further 
manipulation of the iHOMA2 models 
renal excretion threshold variable (R 4 ) 
suggested that alteration of the renal 
threshold by 17% (R 4 = 117%) gave the 
closest approximation for the observed 
mean decrease of 0.9 mmol/L in glucose. 

CONCLUSIONS— 1HOMA2 is an ex- 
tension of HOMA2 that enables the in- 
teractive modeling of physiology and 
treatment effects. It is now possible to 
use more of the potential of the structural 
basis of the model through the 24 in- 
teractive variables. This circumvents a 
drawback of HOMA2, namely, that this 
model does not allow individual compart- 
mental or biological variables to be spec- 
ified even when these are known. 
iHOMA2 can be used in default mode 
(for comparison with all published data 
using HOMA and HOMA2); in analytic 
mode, which allows %B and %S to be cal- 
culated from fasting insulin and glucose; 
or in predictive mode that allows fasting 
insulin and glucose to be calculated from 
%B and %S. 

The iHOMA2 model presents all those 
variables that have some identifiability, so, 
for example, a renal threshold or glomer- 
ular filtration rate (slope of the renal excre- 
tion curve) can be entered to give greater 
specificity to the model. A group that has 
validated data on the observed glomerular 
filtration rate could then validate a specific 
change to the model variable. 

We have demonstrated iHOMA2 use 
in predicting the drug effect of existing 
compounds with a known method of 
action, as in the example of pioglitazone. 
It was shown that of the three hypotheses 
of the site of insulin sensitivity change in 



the basal state, the hypothesis that the 
effect was mainly hepatic was the most 
robust. This is in line with the published 
literature (21,22) and accords with clamp 
data (23) that suggest that a substantial 
effect is mediated by hepatic gluconeo- 
genesis, although the model, with multi- 
ple variables that can be manipulated, 
cannot be regarded as yielding a definitive 
answer to the question. As with all mod- 
els, the conclusion should be taken to be 
illustrative rather than definitive, which is 
further exemplified by the inability of the 
Bland-Altman statistics to discriminate 
between the hypotheses. 

We have exemplified iHOMA2 used 
in modeling with an SGLT2 inhibitor, 
where the effect on blood glucose is 
independent of (3-cell function and insu- 
lin sensitivity. The model suggested that 
the renal threshold change produced by 
dapagliflozin in our study set was best ap- 
proximated by 17% and not 20% as re- 
ported by others. The iHOMA2 model 
can be used as a method of assessing 
quantitative identifiable parameters — in 
this case, renal threshold — from observed 
data. This has the potential, for example, to 
quantify SGLT2 dose effects or compari- 
sons between different SGLT2 inhibitors. 

The utility of a model depends on its 
identifiability. A model is structurally 
identifiable if its parameters can be 
uniquely estimated by fitting the model 
to observed data. This is related to the 
sensitivity of the process outputs to the 
parameter variations and understanding 
how the parameter relationships might 
interact (24-26). The variables in this 
model all have identifiable functions that 
are visualized in each panel as changes to 
the interactive kinetics. This is quite un- 
like curve-fitting exercises, where the 
mathematics may mimic the response 
only because multiple, but physiologi- 
cally unidentified, terms have been used. 
For example, a complex reaction curve 
may be described by a polynomial, but 
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there would be no clear physiological 
"sense" in knowing the value of a squared 
or cubed constant. Model descriptions 
only give insights where the variables rep- 
resent some aspect of the known physiol- 
ogy. With a progressive understanding of 
human physiology, it will be possible to 
be more specific about the model character- 
istics; e.g., one can set (3-cell functional 
changes to be K m or V max related and insulin 
resistance to be tissue specific. Additional 
modules — such as the incretin axis — could 
be activated and manipulated. The resultant 
model could then be used to give predic- 
tions for combinations of agents, activating 
specific mathematical subsections describ- 
ing human pharmacology. 

Mathematical models, such as the 
Mari model (4) or Bergman minimal 
model (3), are useful for the study of com- 
plex phenomena, as they demonstrate 
how measurements of glucose and insulin 
can be seen as manifestations of a com- 
plex feedback process. iHOMA2 has the 
functionality to display the individual 
components of the glucose and insulin in- 
teraction and resolve the steady-state 
mathematics to yield measures of (3-cell 
function and insulin resistance. Using 
the model interactively, one can predict 
(and display graphically) many patholog- 
ical and physiological scenarios. iHOMA2 
can also be individualized, by adjustment 
of the variables, to specific phenotypes, 
genotypes, populations, or individuals. 

The iHOMA2 model does have the 
disadvantage in that, in common with all 
mathematical models, it has limitations to 
its use and is not intended to provide a 
"complete" metabolic description. Allow- 
ing many variables to be internally mani- 
pulated means that identifiability of a 
solution can be at risk (its identifiability 
is reduced). So, iHOMA2 records a list of 
the 24 variables currently able to be 
changed and their settings. For transpar- 
ency in reporting, we would advise that 
the values of the settings on the model are 
explicitly stated when they differ from the 
default. iHOMA2 is still a homeostatic 
model; i.e., it calculates a steady-state so- 
lution. The HOMA2 model was expanded 
to include glucose infusion kinetics (27). 
Further work is underway in iHOMA2 to 
add dynamic aspects so that glucose tol- 
erance tests or clamp data could be mod- 
eled. A drawback of iHOMA2 is that it 
loses the simple assessment aspect of 
HOMA2 — where there are two outputs 
from two inputs. iHOMA2 is, instead, a 
more user-definable model allowing the in- 
dividual characteristics of organs and 



tissues to be manipulated. This manipulation 
allows for hypothesis testing of treatment ef- 
fects and the modeling of pharmaceutical 
treatments for clinical care. 

In conclusion, the creation of the 
iHOMA2 model addresses some of the 
limitations with our previous invariant 
model. iHOMA2 is an application for ex- 
tending the range of the existing HOMA2 
model and can be used to evaluate thera- 
peutic agents in greater physiological detail 
and can predict their effects on fasting 
glucose and insulin and on (3-cell function 
and insulin sensitivity. iHOMA2 can be used 
for hypothesis testing, for the evaluation of 
pharmaceutical agents, for the estimation of 
effect sizes of therapies, and for evaluation of 
indicative changes to glucose when a combi- 
nation of agents is used. It can be used in its 
default mode to generate HOMA2 results 
directly comparable with those in the litera- 
ture. The iHOMA2 program is freely avail- 
able (http://www.ihoma.co.uk) for academic 
and noncommercial use. 



Acknowledgments — The authors thank the 
Oxford National Institute for Health Research 
(NIHR) Biomedical Research Centre for the 
support of N.R.H. D.R.M. acknowledges sup- 
port from the NIHR, of which he is a Senior 
Investigator. 

The authors thank Takeda PLC and Novo 
Nordisk for unrestricted educational grants 
while this work was undertaken. No other 
potential conflicts of interest relevant to this 
article were reported. 

N.R.H. researched data, wrote the manu- 
script, and coded the iH0MA2 program. J.C.L. 
and D.R.M. provided intellectual input into the 
modeling, wrote parts of the paper and dis- 
cussion, and reviewed and edited the manu- 
script. D.R.M. is the guarantor of this work and, 
as such, had full access to all the data in the 
study and takes responsibility for the integrity 
of the data and the accuracy of the data analysis. 

The authors acknowledge the help of 
Dr. R. Newman, University of Oxford, Oxford, 
U.K., for his assistance is creating the online 
iHOMA2 software. The authors also acknowledge 
the work of all the contributors to the previous 
versions of HOMA (6-8) including but not 
limited to the following: J. Hosker, J. Levy, 
D. Matthews, B. Naylor, S. O'Rahilly, A. Rudenski, 
D. Treacher, and R. Turner, with a special 
mention of R. Morris, who extended the model 
considerably to include the proinsulin and 
renal modules. All were members of the Uni- 
versity of Oxford during the development of 
HOMA. 



References 

1 . DeFronzo RA, Tobin JD, Andres R. Glucose 
clamp technique: a method for quantifying 



insulin secretion and resistance. Am J 
Physiol 1979;237:E214-E223 

2. Harano Y, Ohgaku S, Hidaka H, et al. 
Glucose, insulin and somatostatin infusion 
for the determination of insulin sensitivity. 
J Clin Endocrinol Metab 1977;45:1124- 
1127 

3. Bergman RN, Ider YZ, Bowden CR, 
Cobelli C. Quantitative estimation of in- 
sulin sensitivity. Am J Physiol 1979;236: 
E667-E677 

4. Mari A, Pacini G, Murphy E, Ludvik B, 
Nolan JJ. A model-based method for as- 
sessing insulin sensitivity from the oral 
glucose tolerance test. Diabetes Care 2001; 
24:539-548 

5 . Mari A. Mathematical modeling in glucose 
metabolism and insulin secretion. Curr 
Opin Clin Nutr Metab Care 2002;5:495- 
501 

6. Matthews DR, Hosker JP, Rudenski AS, 
Naylor BA, Treacher DF, Turner RC. Ho- 
meostasis model assessment: insulin resis- 
tance and beta-cell function from fasting 
plasma glucose and insulin concentrations 
in man. Diabetologia 1985;28:412-419 

7. Turner RC, Rudenski AS, Matthews DR, 
Levy JC, O'Rahilly SP, Hosker JP. Appli- 
cation of structural model of glucose- 
insulin relations to assess beta-cell function 
and insulin sensitivity. Horm Metab Res 
Suppl 1990;24:66-71 

8. Levy JC, Matthews DR, Hermans MP. Cor- 
rect homeostasis model assessment (HOMA) 
evaluation uses the computer program. Di- 
abetes Care 1998;21:2191-2192 

9. Wallace TM, Levy JC, Matthews DR. Use 
and abuse of HOMA modeling. Diabetes 
Care 2004;27:1487-1495 

10. McAuley KA, Mann JI, Chase JG, Lotz TF, 
Shaw GM. Point: HOMA — satisfactory for 
the time being: HOMA: the best bet for the 
simple determination of insulin sensitivity, 
until something better comes along. Diabetes 
Care 2007;30:2411-2413 

11. Hockaday D, Sayyad M, Yajnik C. Coun- 
terpoint: appreciating homeostasis model 
assessment: more useful earlier rather 
than later. Diabetes Care 2007 ;30 : 24 1 4- 
2418 

12. Hosker JP, Burnett MA, Davies EG, Harris 
EA, Turner RC. Sulphonylurea therapy 
doubles B-cell response to glucose in type 
2 diabetic patients. Diabetologia 1985;28: 
809-814 

13. Degn KB, Juhl CB, Sturis J, et al. One 
week's treatment with the long-acting 
glucagon-like peptide 1 derivative liraglutide 
(NN22 1 1) markedly improves 24-h glycemia 
and alpha- and beta-cell function and re- 
duces endogenous glucose release in pa- 
tients with type 2 diabetes. Diabetes 2004; 
53:1187-1194 

14. Ehrenkranz JR, Lewis NG, Kahn CR, Roth 
J. Phlorizin: a review. Diabetes Metab Res 
Rev 2005;21:31-38 

15. Altman DG. Pratical Statistics for Medical 
Research. London, Chapman & Hall, 1991 



care . diabetesj ournals . org 



Diabetes Care, volume 36, August 2013 



2329 



iHOMA2: an interactive HOMA model 

1 6 . Charbonnel BH, Matthews DR, Schernthaner 
G, Hanefeld M, Brunetti P; QUARTET Study 
Group. A long-term comparison of pio- 
glitazone and gliclazide in patients with 
Type 2 diabetes mellitus: a randomized, 
double -blind, parallel-group comparison 
trial. Diabet Med 2005;22:399-405 

1 7 . Schernthaner G , Matthews DR, Charbonnel 
B, Hanefeld M, Brunetti P; Quartet [cor- 
rected] Study Group. Efficacy and safety of 
pioglitazone versus metformin in patients 
with type 2 diabetes mellitus: a double- 
blind, randomized trial. J Clin Endocrinol 
Metab 2004;89:6068-6076 

18. Komoroski B, Vachharajani N, Feng Y, Li 
L, Kornhauser D, Pfister M. Dapagliflozin, a 
novel, selective SGLT2 inhibitor, improved 
glycemic control over 2 weeks in patients 
with type 2 diabetes mellitus. Clin Phar- 
macol Ther 2009;85:513-519 



19. Cook JT, Page RC, O'Rahilly S, et al. 
Availability of type II diabetic families for 
detection of diabetes susceptibility genes. 
Diabetes 1993;42:1536-1543 

20. Ferrannini E, Ramos SJ, Salsali A, Tang W, 
List JF. Dapagliflozin monotherapy in type 2 
diabetic patients with inadequate glycemic 
control by diet and exercise: a randomized, 
double-blind, placebo-controlled, phase 3 
trial. Diabetes Care 2010;33:2217-2224 

21. Sugiyama Y, Shimura Y, Ikeda H. Effects 
of pioglitazone on hepatic and peripheral 
insulin resistance in Wistar fatty rats. 
Arzneimittelforschung 1 990 ;40 : 43 6-440 

22. Smith U. Pioglitazone: mechanism of ac- 
tion. Int J Clin Pract Suppl 2001:13-18 

23. Gastaldelli A, Miyazaki Y, Mahankali A, 
et al. The effect of pioglitazone on the 
liver: role of adiponectin. Diabetes Care 
2006;29:2275-2281 



24. Cobelli C, Pacini G, Toffolo G, et al. Es- 
timation of insulin sensitivity and glucose 
clearance from minimal model: new in- 
sights from labeled IVGTT. Am J Physiol 
1986;250:E591-E598 

25. Kahn SE, Prigeon RL, McCulloch DK, 
et al. Quantification of the relationship 
between insulin sensitivity and beta-cell 
function in human subjects. Evidence 
for a hyperbolic function. Diabetes 1993; 
42:1663-1672 

26. Jaqaman K, Danuser G. Linking data to 
models: data regression. Nat Rev Mol Cell 
Biol 2006;7:813-819 

27. Hosker JP, Matthews DR, Rudenski AS, 
et al. Continuous infusion of glucose 
with model assessment: measurement of 
insulin resistance and beta-cell func- 
tion in man. Diabetologia 1985;28:401- 
411 



2330 



Diabetes Care, volume 36, August 2013 



care .diabetesj ournals .org 



